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By solving the Bogoliubov - de Gennes equations at zero temperature, we study the effects of a 
one-dimensional optical lattice on the behavior of a superfluid Fermi gas at unitarity. We show that, 
due to the lattice, at low densities the gas becomes highly compressible and the effective mass is 
large, with a consequent significant reduction of the sound velocity. We discuss the role played by the 
lattice in the formation of molecules and the emergence of two-dimensional effects in the equation 
of state. Predictions for the density profiles and for the frequency of the collective oscillations in 
the presence of harmonic trapping are also given. 
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I. INTRODUCTION 

Ultracold gases in optical lattices provide a new fron- 
tier of research where many remarkable phenomena can 
be observed and investigated [l[ . By using Feshbach res- 
onances one can tune the interaction between atoms and 
investigate the BCS-BEC crossover, passing through a 
resonant regime where the scattering length is very large 
and the system exhibits universal properties (unitary 
regime) 2]. First experiments with Fermi superfluids in 
one-dimensional (ID) optical lattices [|| have focused on 
the study of the critical velocity along the crossover and 
revealed that superfluidity is particularly robust at uni- 
tarity. The results are in qualitative agreement with the 
predictions for the Landau critical velocity [H and for the 
behavior of the critical current through a single barrier 

It is well known that a periodic potential favors the 
formation of molecules in a two-component Fermi gas (see 
@ and references therein), providing an effective shift of 
the resonance and bound states in the two-body problem 
even at unitarity. A major problem is to understand the 
consequences of the molecular formation in the superfluid 
phase. Moreover, for large laser intensities, the lattice is 
expected to give rise to 2D effects. 

In order to investigate these properties we use a mean- 
field theorybased on the Bogoliubov - de Gennes (BdG) 
equations [Tj. Although approximate, this approach cap- 
tures basic features along the whole BCS-BEC crossover 
Q, including the formation of molecules and the most 
challenging unitary limit where, for uniform 3D config- 
urations, the predictions arc in reasonably good agree- 
ment with ab initio Monte Carlo simulations js}. The 
BdG equations apply also to situations where the density 
varies over distances of the order of the healing length. 
An important example is given by configurations with 
quantized vorticity [101 ]. Furthermore, BdG equations 



fully account for the modification of the scattering prop- 
erties of fermions induced by the external confinement 
as predicted by Petrov et al. [ll| in the limit of a deep 
periodic lattice. As a consequence, for example, a tight 
ID lattice considerably affects the mean-field superfluid 
transition temperature, where the BCS order parameter 
vanishes [Hj]. Finally, in the case of a deep lattice the 
BdG theory is expected to approach the 2D mean-field 
theory of Ref. [3. 

In this work we study the unitary regime at zero tem- 
perature focusing on the situation in which the lattice 
potential is relatively weak. In this regime, the tight 
binding description is not adequate and a full numerical 
approach based on the BdG equations is called for. By 
solving the BdG equations, we first calculate the equa- 
tion of state, the compressibility and the effective mass 
of the unitary Fermi gas in the lattice. The results are 
then used to obtain interesting predictions for observ- 
able quantities such as the sound velocity, the frequency 
of collective modes, and the density profile of trapped 
gases in typical experimental configurations. 



II. FORMALISM 

At zero temperature, the chemical potential \i of a su- 
perfluid Fermi gas in a lattice is given by the derivative of 
the energy density e = E/V with respect to the average 
(coarse-grained) density n: 



de{n,P) 
dn 



(1) 



where P is the quasi-momentum of the superfluid along 
the lattice [14j. The compressibility k and the effective 
mass m* , are given by the second derivatives of e with 



2 



respect to n and P: 

_!_ d 2 e(n,P) _ dn(n,P) 1 _ 1 d 2 e(n,P) 
<9n 2 9n m* n <9P 2 

(2) 

We calculate these quantities at unitarity for P = 0, i.e. 
for a gas at rest, in the periodic potential 



V ext (z) = sE R sin 2 q B z. 



(3) 



Here s is the laser intensity, E-r = ti 2 q 2 i /2m is the recoil 
energy, = ft/d is the Bragg wave vector, d is the lat- 
tice constant and m is the atom mass. Then the BdG 
equations are given by 



H'(r) A(r) 
A*(r) -H'(r) 



Ui(r) 
Vi(r) 



Ui(r) 
Vi(r) 



(4) 



where H'(r) = — h 2 X7 2 /2m + V ex _t — The order pa- 
rameter A(r) and the chemical potential [i, appearing in 
Eq. are variational parameters determined from the 
self-consistency relation 



A ( r ) = -5 5>*( r )<( r )' 



(5) 



together with the constraint n = (2/V) J |^i(r)| dr, 
enforcing conservation of the average density n. In 
Eq. ([5]) , g is the coupling constant for the contact interac- 
tion and the BdG eigenfunctions obey the normalization 
condition / dr [u*(r)u 3 (r) + v*(r) Vj (r)} = 5 itj . 

For contact interactions, the right hand side of Eq. ([5]) 
is ultraviolet divergent and must be cured by the pseudo- 
potential method. A standard procedure consists of in- 
troducing a cut-off energy Eq in the sums over the BdG 
eigenstates and of replacing the bare coupling constant 
by the s-wave scattering length a s through the rela- 
tion (kpas)- 1 = %TrE ¥ /{gkl) + (2/tt) v / Ec/E f , where 
hp = (S^n) 1 ^ 3 and Ep — h 2 kp / (2m) are the Fermi wave 
vector and energy, respectively, of a uniform noninteract- 
ing Fermi gas with density n. An alternative regulariza- 
tion scheme is the one proposed in Refs. [HI, [l6| . This 
method exploits directly the short distance behavior of 
the single-particle Green's function and is efficient even 
in the presence of a tight confinement. Both procedures 
give identical results provided Eq is large enough. 

In the presence of a supercurrent with wave vector 
Q = P/h moving along the lattice, one can write the 
order parameter in the form A(r) = e t2 ® z A(z), where 
A(z) is a complex function with period d. Therefore, 
from Eq. (JSJ), we see that the eigenfunctions of Eq. 
must have the Bloch form Ui(r) = Ui(z)e l Q z e lk ' T and 
Vi(r) = Vi(z)e~ l Q z ^e* k ' r ', where k z lies in the first Bril- 
louin zone and iii and Vi are periodic in z with period 
d. This transformation reduces Eq. (0| to the BdG equa- 
tions for Uj and Vj as 



H' Q (z) A(z) 
A*(z) ~H'_ Q (z) 




Ui(z) 
Vi(z) 



where 

H' Q (z) = ^-\k 2 x + k 2 + (-id, + Q + k z f] +V ext (z)-fi . 
^ 2m L " J 

(7) 

From now on, the label i represents the wave vector k as 
well as the band index. 

From the solution of the BdG equations, we can di- 
rectly calculate /i and k [see Eqs. (flj and while m* 
is obtained from the energy density e(n,P) 



= dr 



, (6) 



^2( M - e,)l^)| 2 + &*(z)tii(z)Z*(z) 

i i 

(8) 

Note that in the rhs of Eq. ([8]) both contributions are 
separately divergent but the sum is finite, as one can 
easily check for the uniform case (s = 0). 



III. RESULTS 

A first important remark concerns the low density limit 
of n and m* whose values are found in perfect agreement 
with the results of the exact solution of the two-body 
problem Q. This proves that the BdG theory correctly 
accounts for the deep modifications of the atomic scatter- 
ing properties induced by the external confinement which 
gives rise to bound molecules even at unitarity. 

Our results for the density dependence of k _1 and m* 
at unitarity are shown in Figs.[IJa) and[TJb). They are 
plotted for s = 1, 2.5, and 5, as functions of E-p/En = 
(kpd/-K) 2 . The solid, dashed and dashed-dotted lines cor- 
respond to the values obtained with the regularization 
scheme of Ref. [l6[ and choosing the cut-off energy large 
enough to ensure convergence within 1-2% [l7l |. In or- 
der to appreciate the role of the cut-off energy, we also 
plot the results obtained using the standard regulariza- 
tion procedure with Eq = IOOEf (points). The results 
agree within 10-15%, the accuracy being worse in the 
small density regime Ep E^. 

In the absence of the lattice (s = 0) the results take a 
universal behavior at unitarity: the only relevant length 
being the interparticle distance fixed by kp . Due to trans- 
national invariance, one can write e(n,P) — e(n,0) + 
nP 2 /2m so that m* = m. Furthermore, the energy den- 
sity at P — can be written as e(n, 0) = (1 + f3)e°(n, 0), 
where e°(n, 0) = (3/5)71^ is the ideal Fermi gas value 
and [3 is the dimensionless universal parameter account- 
ing for the interactions in uniform gas. The BdG equa- 
tions predict f3 ~ —0.41 to be compared with the Monte 
Carlo result j3 ~ -0.58 0. 



A. Unitary Fermi gas in a lattice 

1. Equation of State 

New features appear when s increases. Let us first 
discuss the behavior of the equation of state. At small 
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expanding in the small parameter sE^/Ep, we find [l£ 
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FIG. 1: (Color online) Inverse compressibility kT 1 , and ef- 
fective mass ra* of the unitary Fermi gas for s = 1 (red), 
2.5 (blue), and 5 (green). Lines are convergent results ob- 
tained by the regularization scheme of Ref. fig. Points are 
obtained by means of the standard regularization scheme with 
Ec = 100-Ef . Asymptotic expressions @ and (|10p are shown 
by the dotted lines. Open circles in panel (b) show m* ob- 
tained in Ref. The s = 1 results for m* are also shown in 
the inset in the linear scale. 



densities (Ep/En <c 1) we find that the lattice causes a 
linear density dependence of the chemical potential [see 
inset of Fig. 2] and a strong suppression of k~ 1 with re- 
spect to the uniform value [see Fig. Ufa)]. These results 
are consistent with the tendency of the lattice to favor 
the formation of molecules. The size of these molecules 
is fixed by the values of s and d so that, in the limit 
of a dilute gas (kpd -C 1), the interparticle distance 
can be larger than the molecular size. In this limit, 
the BdG equation describes the formation of a molec- 
ular gas which gives rise to Bose-Einstein condensation; 
the equation of state is given by the Bogoliubov theory 
and re -1 is expected to be linear in the density, yielding 
k~ 1 /k~ 1 (s = 0) cx rt 1 / 3 — » 0. For s ^> 1 the chemical po- 
tential remains almost linear in density even at relatively 
large densities, due to 2D effects caused by the bandgap 
in the longitudinal motion (see discussion below). At 
even higher densities, one eventually recovers the behav- 
ior of a uniform gas since the lattice only provides a small 
perturbation. By using an hydrodynamic approach and 
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(sE R /E F r (9) 
This is shown by dotted lines in Fig. []Ja). 
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FIG. 2: (Color online) Comparison between the unitary 
Fermi gas (solid lines) and the ideal Fermi gas (dashed lines) 
in the same optical lattice with s = 2.5. The region between 
the two vertical lines represents the range where the chemical 
potential of the ideal gas lies in the bandgap above the lowest 
Bloch band. In panel (b), the square at m*/m = 1.20 shows 
the value of the effective mass of a single atom with mass m 
in the same lattice, while the triangle at m* /m = 5.53 shows 
that of a point-like molecule with mass 2m. 



2. Effective Mass 

The formation of molecules has important conse- 
quences also for to*. Due to this effect, at low densities, 
the enhancement of to* caused by lattice is much more 
dramatic in the unitary Fermi gas compared to the ideal 
Fermi gas (or, equivalently to m* in the BCS limit) and 
to the gas of bosons with the same mass to (see Fig. 4 
in Ref. [lj|). We also note that, for a given s, the value 
of to* at E F lies between the value of to* for a sin- 
gle atom of mass m in the same lattice [square in the 
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vertical axis of Figs. [2tb)] and the value of to* calcu- 
lated for a point-like molecule of mass 2m (triangle) . As 
Ey/Eh increases, the effective mass exhibits a maximum 
at E-p/Es. ~ 1, then it decreases towards the bare mass, 
to* = to. The same hydrodynamic argument used for 
k _1 also explains the behavior of the effective mass for 
small sE^/E-p: 



m 
m 



-0 



{sE R /E F f 



> (io) 



showing that the effect of the lattice is stronger for to* 
than for kT x . It is worth comparing the results with the 
case of bosonic atoms, where to* decreases monotonically 
with increasing density since the interaction broadens the 
condensate wave function and favors the tunneling [l9j . 



3. Comparison with the Ideal Fermi Gas 
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FIG. 3: (Color online) Sound velocity c s in units of the 
sound velocity d 0) = [(2/3)(l + f3)E F /m] 1/2 for the uniform 
system. As in Fig. ffj red, blue and green lines correspond to 
s = 1,2.5, and 5, respectively. 



The occurrence of a maximum in the curves for both 
kT x and to* can be interpreted as an effect of the energy 
gap in the longitudinal motion, which opens at q = qb- 
An instructive comparison can be made with an ideal 
Fermi gas in the same lattice, where the effects of the 
bandgap are more evident due to the sharper Fermi sur- 
face. In Fig. [J] one sees that the ideal gas curves have 
two cusps. They occur precisely when ji coincides with 
the top of the lowest band and the bottom of the first 
excited band at q = q-g, respectively. In between, as 
Ep/Eji increases, only the transverse modes are avail- 
able and the system behaves effectively like a 2D system, 
the longitudinal degree of freedom simply giving a con- 
stant contribution to kT 1 and to*. Consequently, n~ l 
and to* are proportional to the average density in this 
region (see Eqs. (10) and (11) in Ref. [2fj]) and one finds 
the power laws k~ 1 /k^ 1 (s = 0) oc Ep and to* oc E f ^ 2 . 
The interval where this 2D behavior takes place becomes 
wider as s increases, since the bandgap increases with s. 
By comparing the ideal and unitary Fermi gases, we see 
that the interaction significantly smears the effect of the 
gap as a result of a much broader Fermi surface. Espe- 
cially the maximum values of to* are drastically reduced. 
Note also that the molecular- like pair correlations, which 
are responsible for the low density behavior of k _1 and 
to*, are absent in the ideal Fermi gas. 

4- Sound Velocity 

Our results for k" 1 and to* can be used to calculate the 
sound velocity c s = y k _1 /to*, which is given in Fig. [3] 
One can see a significant reduction of c s compared to the 
uniform system. Except for very low densities, the sound 
velocity varies rather smoothly with the density. This 
makes it possible to provide an estimate the change of 
c s induced by the lattice even for harmonically trapped 
gases, where the coarse-grained density is not uniform. 



Notice that the propagation of sound is a direct conse- 
quence of superfluidity and, in the presence of tight lat- 
tices, can be regarded as a Josephson effect, where the 
gas tunnels in a coherent way through the barriers pro- 
duced by the lattice. 



B. Unitary Fermi gas in a lattice + trap 

1. Density Profile 

We can also provide useful predictions for the density 
profile and for the collective motion of a trapped gas, 
when a harmonic confinement is added to the periodic 
potential. The coarse-grained density profile, n(z), is 
easily calculated using the local density approximation 
(LDA) for /i. Figure S] clearly shows that, for s = 5, the 
profile takes the form of an inverted parabola, reflecting 
the linear density dependence of the chemical potential 
(see inset). In this calculation, we set lu± = lu z , where 
u>± and lu z are the transverse and longitudinal trapping 
frequencies, hhJ z / 'En = 0.01, and the number of particles 
N = 10 6 ; these parameters are close to the experimental 
ones in Ref. 

The accurate measurement of the density profile in the 
presence of harmonic trapping can actually give direct 
information on the compressibility of the gas. In fact, 
within the range of validity of LDA, the two quantities 
obey a simple relationship which, for isotropic traps, is 
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FIG. 4: (Color online) Coarse-grained density profiles of a 
trapped gas, n(r± — 0, z) for s = and 5 in units of the 
central density n(0) = 0.0869^1 calculated for s — (this 
local density corresponds to Ef/Er, = 1.88). The quantity 
is the axial Thomas-Fermi radius for s = 0. The inset 
shows the density dependence of the chemical potential. The 
parameters of the trap are given in the text. 



function of n as one moves from the center to the border 
of the atomic cloud. For clouds whose maximum den- 
sity is such that E-p/En < 1.5 the minimum value cor- 
responds to the low density Ev/En — » limit, while the 
maximum corresponds to the absolute maximum of the 
curves in Fig. [ljb). In this way we obtain the estimate 
0.89 < lu d /lu z < 0.96 for s = 1, 0.59 < uj d /lo z < 0.66 for 
s = 2.5, and 0.26 < uj d /lu z < 0.28 for s = 5. 

Finally, the effects of the lattice on the equation of 
states can be observed also by studying the compression 
modes. For example, taking the cigar shape geometry 
lo z <C the frequency of the radial breathing mode 
is unaffected by m* . In the regime where the chemical 
potential is linear in the density (see inset of Fig. [4|) the 
collective frequency approaches the value uj = 2<x>j_, while 
its value is y/lQ/Zu>±_ in the absence of the lattice. The 
transition between the two different regimes is similar to 
the one recently investigated by varying the scattering 
length on the BEC side of the resonance in the absence 
of periodic potentials (see, e.g., 0). 



IV. CONCLUSION 



2. Collective Modes 

The collective modes can be studied by solving the 
hydrodynamic equations, where the periodic potential is 
included through its effects on the equation of state n(n) 
and on the effective mass which determines the current 
in the longitudinal direction 2l|. While the frequency 



of the dipole (center-of-mass) oscillation along the trans- 
verse direction is not affected by the lattice, a useful es- 
timate for the dipole frequency in the z-direction can 
be obtained using a sum rule approach based on the 
calculation of the energy weighted and inverse energy 
weighted moments of the dipole strength within the hy- 
drodynamic theory 20]. This approach yields the result 
u>D = uj z (m/ to*) 1 / 2 where 



1 



1 

N 



dr 



n(r) 



i(r)] 



(12) 



Even without calculating n(r), we can easily estimate 
lower and upper bounds for the frequency by replacing 
to* with the maximum and minimum values of m* as a 



In this work we have shown that the inclusion of a 
ID optical lattice, by favoring the formation of molecu- 
lar configurations and by inducing a band structure in 
the quasiparticle spectrum, has profound consequences 
on the thermodynamic quantities, the density profile, and 
the collective oscillations of the unitary Fermi gas. Our 
calculations can be naturally extended to the case of fi- 
nite quasi-momenta, where energetic and dynamic insta- 
bilities are expected to impose some limits to the super- 
fluid motion. Further investigations of the 2D nature of 
the many-body system achievable with large laser inten- 
sities are also in progress. 
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